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ABSTRACT 

A very important issue concerning protostellar jets is the mechanism behind their formation. Ob¬ 
taining information on the region at the base of a jet can shed light into the subject and some years ago 
this has been done through a search for a rotational signature at the jet line spectrum. The existence 
of such signatures, however, remains controversial. In order to contribute to the clarification of this 
issue, in this paper we show that the Principal Component Analysis (PCA) can potentially help to 
distinguish between rotation and precession effects in protostellar jet images. This method reduces 
the dimensions of the data, facilitating the efficient extraction of information from large datasets as 
those arising from Integral Field Spectroscopy. The PCA transforms the system of correlated coor¬ 
dinates into a system of uncorrelated coordinates, the eigenvectors, ordered by principal components 
of decreasing variance. The projection of the data on these coordinates produces images called to¬ 
mograms, while eigenvectors can be displayed as eigenspectra. The combined analysis of both can 
allow the identification of patterns correlated to a particular physical property that would otherwise 
remain hidden, and can help separating in the data the effect of physically uncorrelated phenomena. 

These are for example, rotation and precession in the kinematics of a stellar jet. In order to show the 
potential of the PCA analysis, we apply it to synthetic spectro-imaging datacubes generated as an 
output of numerical simulations of protostellar jets. In this way we generate a benchmark to which a 
PCA diagnostics of real observations can be confronted. Using the computed emission line profiles for 
[O I]A6300 and [S II]A6716, we recover and analyze the effects of rotation and precession in tomograms 
generated by PCA. We show that different combinations of the eigenvectors can be used to enhance 
and to identify the rotation features present in the data. Our results indicate that the PCA can be 
useful for disentangling rotation from precession in jets with an inclination of the jet with respect to 
the plane of the sky as high as 45°. We have been able to recover the initially imposed rotation jet 
profile for models at moderate inclination angle {(j) < 15°) and without precession. 

Subject headings: ISM: Herbig-Haro objects — ISM: jets and outflows — ISM: kinematics and dy¬ 
namics — methods: data analysis 


1. INTRODUCTION 

Recent observations of protostellar jets show, in some 
cases, a systematic side-to-side (with respect to the jet 
axis) shift in radial v elocities in se v eral e missi on lines in 
the optical a nd UV. Davis et al.| ([2000) and Bacciotti 
et al.| (|2002|) reported the first results of this kind, in 
a molecular and atom ic jet, respectively. In particular, 
Bacciotti et al. (2002) presented an analysis performed 
on high resolution Hubble Space Telescope (HST) long- 
slit spectra taken at different positions, both along and 
across the DG Tau jet, in the region close to its driving 
source. They found a systematic side-to-side difference 
in the radial velocity of the jet with respect to its axis, 
which they interpreted as evidence of jet rotation, what 
was later found to be in agreement with the sense of 
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rotation of the circumstellar disk of DG Tau (Testi et al. 


2002 ). 


This observational finding has been followed by other 
studies which tried to look for the same kind of pat¬ 
tern in other jet sources. Goffey et al. (2004|), using 
the Space Telescope Imaging Spectrograph (STlS) of the 
HST found similar results for the jets in the T Tauri 
stars TH 28, RW Aur and LkRa 321. Gonsistently with 
the predictions of the magnet o cent rifugally driven mech¬ 


anism for jet launching (e.g., [Blandford & Payne [1982 
Ferreira||1997), the bipolar collimated outflows in TH 2^ 
and RW Aur showed the same signs of velocity asymme¬ 
try in both red- and blue-shift lobes, implying the same 
sense of rotation in both. However, the sense of the jet 


of the disk rotation (Cabrit et al. 

2006|). Furthermore, 

recent observations by Gottey et aJ 

. (|2U12) showed a ra- 


corresponding to a rotation in the opposite sense of that 
determined by the optical observation. In addition, no 
shift in velocity was detected in the observation of RW 


Aur with t he same instrument six months later (Coffey 

et al.||2012| ). 


T he CW Tau and HH 30 jets were observed in the near 
ult raviolet (NUV) and optical wavelengths, using STIS 
by Coffey et al. (|2007), who also presented new observa- 
tions of the DG T'au and TH 28 jets. A systematic radial 
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velocity shift pattern was observed in the CW Tan jet 
spectra in the same range of velocities, of ^ 10 — 20 km 
s~^ for the optical shifts, with slightly smaller shifts in 
the NUV (^5 — 10 km s“^), while no significant shift 
in radia l velocity was obs erved for the HH 30 jet. For 
this jet, |Pety et al. (2006|) did not observe also any ro- 
tati on signature at m iillimeter wavelengths. As noted by, 
e.g.,|Cai et al.|(|2008|), the lack of any rotational evidence 
for HH 30 jet is unexpected since the rotation signature 
should be maximum for jets moving in the plan e of the 
sky as is the case of HH 30 (inclination angle < 1 


The indication of rotation in DG Tan, TH 28 and 
LkH(al28 is also consistent with a wind launched at the 
innermost part of the accretion disk. It also reproduces 
the onion like structure (in radial velocit ies) expected for 
a wind launched by an accretion disk. [Anderson et al 


(2003) found for the DG Tan jet, using Bacciotti et al 
(2002) data, a launching radius of vq ^ 0.3 — 4 AU f rom 
the YSO, thus excluding the X-wind model (e.g., Gai 


Mundt 


These optical and NUV observations of a sample of jets 
have been followed by near-infrare d (NIR) long sl i t spec - 
troscopy of HH 212 carried out by Gorreia et al. (2009). 
HH 212 is one of the bona fide examples of jet symmetry 
(between both lobes) and th e first HH jet for w hich evi¬ 
dence ofrotatmm^sfomd (Davis et al. 12000). In their 

4l( 


work, 
feet o: 


Gorreia et al. 


([2009) argued that the combined ef- 
rotation and precession of the jet axis would be 

responsible for the observed pattern, and that rotation __ 

alone would not be able to account for the data. The HH sured (e.g., Th28, see Goffey, Bacciotti & Podio 2008) 


et al. 2008) as a jet launching mechanism. The obser 
vational determination of the launching region is, unar- 
guably, a key piece to advance further and to gather all 
the information in a consistent and complete model. 

However, the interpretation of side-to-side velocity 
shifts as evidence of rotation is still controversial. Side to 
side asymmetries in the velocity along the main jet axis 
of order of 10% of the observed jet speed would produce 
an effect on the radial velocity similar to that observed 
in the works cited above. In this sense, the lack of veloc¬ 
ity shifts in HH 30 and the presence of clear side-to-side 
asymmetries in the electron and hydrogen densities of 
some of the jets where the radial velocity shifts are mea- 


212 jet has also been observed by Godella et al. (2007), 
who did not find evidence of velocity gradients compati- 
ble with rotation in SiO observations. 

Along with these efforts, attempts have been made to 
observe rotation of t he jet axis through ground based, 
world class telescopes. [Goffey et al.|(|2011|) used the Gem¬ 
ini South Telescope to observe HH 34, HH 111 and HH 
212. They found evidences for rotation in some knots of 
HH 111 and 212. For HH 111, they noted that the sense 
of rotation obtained for the observed knot is in opposite 
sense to the disk rotation, while for HH 212 they agree 
(in bc )th receding and app roaching lobes, and the disk it¬ 
self). Goffey et al. (2015) looked for radial velocity shifts 
in the KY 'Fau system, by using Gemini NIFIS+ALTAIR 
observations of the jet, combined with radio observations 
of the disk using the Plateau de Bure. Although a Kep- 
lerian rotation pattern for the disk was clearly obtained, 
the value of the radial velocity shifts remained below the 
3cr detection limit. However, the authors considered the 
obtained radial velocity values as upper limits and as¬ 
suming steady state, constra i ned th e launching region to 


are consistent with the possibility that at least part of 
the observed radial velocity shifts originates from side- 
to-side asymmetries in the velocity of the jet along the 
d irection of the jet propa gation. 


Gerqueira et al. (2006) discussed the possibility that 
the precession of the jet axis could also give rise to the 
same radial velocity asymmetries ob s erved on jet spec¬ 
tra by, for example, ~ 
investigated the com' 


Bacciotti et al. (2002). They also 
Dined ettect of precession and rota¬ 
tion, which seem to operate together in the D G Tau jet. 


be below 0.45 AU. In Fendt (2011) a comprehensive list 
and a detailed discussion of HH objects for which there 
is evidence for jet rotation are presented. 

If the rotation interpretation is correct, these stud¬ 
ies open an important new window to investigate the 
physics behind the jet production and launching mecha¬ 
nisms. The measurement of jet rotation has implications 
in the estimate of the angular momentum flux lost by 
the system, on the determination of the footpoints of 
the jet launching region, and could help to discriminate 
among differe nt jet production models (see, for instance, 
Pudritz 2004[). The comparison between the expected 
values for the radial velocity asymmetry, from magneto¬ 
hydrodynamics (MHD) disk wind models , and t he values 
obtained f or DG Tau by [ Baccio tti et al.[ ([2002[) was pre¬ 
sented by [Pesenti et al.[ ([2UU4[). They found that both 
classes of MHD disk wind models, the so-called cold and 


warm solutions (see Gasse & Ferreira[[2000a|b), are able 
to adjust the observed trend for DG 4'au transverse radial 
velocity, but only the warm solution is able to reproduce 
the velocity shifts. 


While it is c lear that the models presented by [Gerqueira 
et al. (2006) do not apply for je ts with small (5, 5° ) pre¬ 
cession angles, as mentioned by[Goffey et al.[([2007[), it is 
also evident that at least for those systems that show ev¬ 
idence for precession (e.g., DG Tau), its presence should 
be considered properly. 

We propose here a method that could be used to distin¬ 
guish between rotation and precession and in general any 
asymmetry in the velocity components of the jet when 
interpreting a given spectra from a Fabry-Perot or an 
Integral Field Unit (IFU)-hke data. For this purpose, 
we perform a Principal Gomponent Analysis (PGA) of 
synthetic spectra generated from the nu merical simula¬ 
tions performed in [Gerqueira et al.[ ([2006[) , intended here 
as idealized reproductions of the observations. Following 
this approach, we construct 3D synthetic datacubes in 
two relevant emission lines ([S II]A6716 and [O I]A630(Q 
and then apply the PGA technique to split the spectra 
into its components ranked by the variance. The main 
goal of this work is to show how the PGA can help dis¬ 
entangling the different mechanisms that cause asymme¬ 
tries at the jet lines for rotating and processing jets. 

The paper is organized as follows. In Section [^ we 
discuss the PGA technique as applied in our numerical 
simulations, which are briefly reviewed in Section [^ In 
Section [^ we discuss tomograms and eigenspectra and 

® We have actually applied the PC A technique also to datacubes 
from [O l]A6363 and [S Il]A6731 emission lines. We found, how¬ 
ever, that the results for the two oxygen lines are the same, and 
that the results for the two sulfur lines are very similar. There¬ 
fore, throughout the paper we will discuss the results for only one 
emission line for each atom. 
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their physical interpretation. In Section the results 
based on the reconstructed process of the treated dat- 
acube are shown and in Section [6] the main conclusions 
are presented. In the Appendix [A] w e present Tables 
for eigenvalues and in the Appendix |B| we show how the 
analysis change when we take into account noisy data 
and how we can extract the desired information through 
the use of the enhancement factor in PCA. 


2. THE PCA TECHNIQUE APPLIED TO A DATACUBE 

The Principal Component Analysis (PCA) is a method 
to extract information from large, multidimensional 
datasets, to identify peculiar patterns in the data that 
would otherwise remain hidden or mutually combined. 
The isolated patterns can in many cases be associated 
to physical properties that would remain undetected in 
traditional spectro-imaging diagnostics. This is done 
through the construction of a new system of “natural” 
uncorrelated coordinates ordered by decreasing variance 
with respect to the average image/spectrum. This new 
system describes the dataset in a more efficient way, as 
the projection of the data on each coordinate isolates one 
relevant information. 

In the presentation of the PCA technique , we will use 
the sa me fo rmalism of Steiner et a l. (2009), Ricci et al. 
(2011) and Menezes et al. ( 1 - 2014 . These authors de¬ 
veloped a PCA-based tool to investigate the properties 
of datacubes obtained from observations made with the 
Gemini Multi-Object Spectrograph in its Integral Field 
Unit mode (GMOS-IFU) for active galaxies (NGC 4736, 
NGG 7097 and NGC 3115, respectively). First, they use 
PCA, together with other cleaning techniq ues, to remove 
instrume ntal spatial fingerprints (see also |Cerqueira et 
[aL]|2Q15| who applied the same cleaning technique to 
uMOS-IFU observations of the HH 111 jet). As we shall 
see below, in the PCA technique a given datacube is 
firstly reduced to a particular bi-dimensional array from 
which one calculates the covariance matrix. Eigenvectors 
and eigenvalues are then determined for this matrix, and 
ordered by decreasing variance. The eigenvectors consti¬ 
tute the new orthogonal basis, and the data can be pro¬ 
jected on this basis to form the tomograms. Each tomo¬ 
gram will highlight a particular pattern in the data. The 
tomograms^ that are bi-dimensional maps which describe 
the data in an orthogonal basis of uncorrelated coordi¬ 
nates, can then be interpreted and can reveal physical 
properties associated with the obtained pattern. Tomo¬ 
grams can be combined linearly to reconstruct cleaned 
2D images in the normal space in which that particular 
property is now evident. 


2.1. Preparing the data for the PCA analysis 

In the more general case, we simulate a jet that can 
precess and rotate (see Eigurej^. The axis of precession 
makes an angle (f with respect to the plane of the sky, 
assumed here to be the (x,y) plane. The axis of the 
rotating/precessing jet lies along a cone that makes an 
angle 0 with respect to the axis of precession. A non 
processing, but rotating jet has 0 = 0. We compute as 
an output of the numerical simulation the flux of each 
emission line arranged in a tridimensional datacube with 
two spatial and one spectral dimensions. Each element of 
the datacube collects the flux produced in a square region 
of 5.74 AU (the spatial resolution of the simulation) in 


(x,y), in the radial velocity interval from -400 km s ^ 
to -j-lOO km s“^ (corresponding to a spectral width of 
- 10 A for [O I]A6300 and - 11 A for [S II]A6716), 


with a spectral sampling of 10 km s“^. The datacube 
will have n = ja x n = 128x128 pixels in space and m = 
50 pixels in wavelength. Our data is equivalent to those 
generated by an lEU: from the simulation, we can obtain 
velocity channel maps around a given emission line for a 
fixed wavelengt h. By varying the wav elength, we obtain 
a datacube (see Cerq^eira et al.||2006 ). 

Eor each value of 
with i = 1,... ,/i; j 


rque 


we can define an intensity, (IijA)o 5 
1,..., z^; A = 1,..., m. The “0” 
index ref ers to the original int ensity extracted from the 
data (see Ricci et al.||2011 and Steiner et al.||2009), i.e., 
the original intensity at each velocity channel map. Since 
/r and n are the total number of pixels in the x and 
y spatial directions we can define the average spectrum 
of the datacube as the intensity averaged over the two 
spatial dimensions: 


= ( 1 ) 

i=l j=l 

where n = jix n. This mean emission is used to redefine 
the datacube as: 


AjfA — (AjfA)o Q\ 


( 2 ) 


which is the excess (or deficit) of flux in comparison to 
the average flux at each pixel at a given wavelength. 

Since we are looking for patterns of variation in the 
data, the subtraction of the averag e spectra is an impor¬ 
tant procedure. As pointed out by Steiner et al.|(2009), 
this also eliminates components of emission in the spec¬ 
tra that have null spatial variance at a given wavelength, 
as is the case of the sky emission that could be eventually 
constant in the observed field in real (i.e., not synthetic) 
datacubes. 

We can now transform the datacube in an intensity 
matrix^ Igx, of n rows (the spatial pixels) and m columns 
(the spectral pixelslq The spatial pixel in the intensity 
matrix (P) is relatedto the datacube spatial pixels {i and 
j) through the relation 




(3) 


This univocal relation is used to i) build up the inten¬ 
sity matrix (I/3a) mentioned before and ii) to recover 
the d atac ube intensity for a given projection (see Sec¬ 
tion ^2.2| below), which can contain the whole data or 
just a fraction of it. In the simulations that we are going 
to present in this paper, ji = v = 128 and m = 50 (and 
so n = 16 384). 

The matrix Igx is then ready to be processed using 
the PCA technique. In the n ext section we discus s the 
PCA pipeline as described in Steiner et al. (2009). We 
present here the same approach but for simulated dat¬ 
acubes. Our aim, as already emphasized, is to use con¬ 
trolled inputs, including different physical mechanisms 

^ The wavelength is transformed in radial velocity with respect 
to the rest fr ame of the star in our velocity channel maps. In |Steiner| 


|et al. ( |2009| > , m refers to the wavelength, which is equivalent. 

° 'I he spatial pixel s represent the objec ts while the spectral pixels 
their properties (see [Steiner et al.|2009| >. 
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like precession and rotation, to see their imprints in the 
tomograms. Here, we want to test the potential of PCA 
in disentangling the effects of rotation and precession in 
real data of protostellar jets. 


2.2. An orthogonal basis of uncorrelated coordinates 

The main goal of the PCA technique is to represent 
the data in a new set of mutually orthogonal basis formed 
finding the eigenvectors and eigenvalues of the covariance 
matrix of the modified dataset Igx, defined as: 


Ccov 




i/3A 


n — 1 


(4) 


(see Steiner et al.||2009) where [I/3 a]^ is the transpose of 
the intensity matrix, Igx. The covariance matrix has m 
eigenvalues, representing the variance of the data for the 
associated eigenvector (where k is the order of the 
eigenvector), that can be ranked in the new basis with 
decreasing value of the variance, i.e. from the largest 
(Al) to the smallest (A^) variance, as just the ones cor¬ 
responding to the largest variances will contain relevant 
information. The variances can be represented in a nor¬ 
malized way in terms of percentages, and they are the 
diagonal elements of the covariance matrix; the sum of 
all variances must return X^^i Aj^ = 1 (or 100%). 

The eigenvectors represent a new orthogonal basis in 
which the data can be described. These eigenvectors 
are used to build the characteristic matrix^ Ea/c, which 
has m eigenvectors in its columns, sorted by decreasing 
variance. In this new orthogonal basis, the data can be 
decomposed as: 


T/3/c = I/3A • Ea/c , (5) 

where ’Tgk is the intensity matrix expressed in the new 
basis. In this basis, the covariance matrix Dcov? defined 
as: 


D, 




L/3/c 


n — 1 


(6) 


is diagonal by construction (orthogonality property), 
which means that it has null covariance between differ¬ 
ent coordinates and its diagonal elements are the eigen¬ 
values. Since the basis is orthogonal, the projection Tgk 
on the space of eigenvectors represents a pattern that is 
uncorrelated for each k, and one can refer to uncorrelated 
physical properties that can be separated and recognized. 

A way to visualize the projection Tgk of the data onto 
the new basis is to produce from it 2D images called to¬ 
mograms. The tomogram is the quantity Tij^k obtained 
by Tgk by selecting a given k = k' and by unpacking p 
into the two indexes by inverting equation It is 
important to emphasize that the tomograms are images 
of the 3D data in the new system of coordinate. As in 
|Steiner et al.|(|2009|), the tomogram for the /c—eige^ector 
IS an image that arises when we use equations and 
for a given /c, to assign a value for the intensity for 
each pixel. The generated image will have the same spa¬ 
tial dimension of the original one. Also, it will be inde¬ 
pendent of A (or the radial velocity), since this “index” 
is summed up in equation Each tomogram corre¬ 
sponds to an eigenspectrum, which are essentially 


TABLE 1 
Simulated models 


Model 

(yr) 

(yr) 

(°) 

Rotation? 

Notes 

Ml 

8 

- 

- 

no 

Reference model 

M2 

8 

8 

5 

no 

Precessing 

M3 

8 

- 

- 

yes 

Rotating 

M4 

8 

8 

5 

yes 

Processing/rotating 


^Tp is the pulsation period of the jet variability, 
is the period of precession of the jet axis, 
is the angle the jet makes with the precession axis (the preces¬ 
sion angle). 


the components of the eigenvector plotted against the 
wavelength (or radial velocity). The tomogram should 
be analyzed together with its associated eigenspectrum, 
which can help to identify the physical process behind 
the observed patterns. In the Section]^ we present and 
discuss several tomograms and eigenspectra. 

The datacube in the real space {iA^ A) can be recon¬ 
structed by inverting the equation 


Ii,i,A = I/3A = T/3 /c • [Ea/c]^ . (7) 

where we have used also the equation (§. 

The characteristic matrix, Ea/c, can be manipulated to 
keep, for example, just the eigenvectors of higher rele¬ 
vance, suppressing the remaining ones. The usefulness 
of such a reconstruction is the possibility to obtain a 
datacube in which some physical properties can be em¬ 
phasized, and some undesirable features presented in the 
data, like in strumental fingerprints or even noise can 
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be removed d Steiner e t ah |2009 Menezes et al. 2014 
I2OI5I). This procedure is well docu- 
r et al. (2009), but we will rewrite their 


mented in Steiner 
equation (8) here for the' sake' of clarity: 


La(< ^max ) = T/3fe(< k max) [Ea/c( ^ ^max)] 


( 8 ) 


In this equation, /Cmax is the limit of the relevance of the 
eigenvector that the user wants to keep in both, the char¬ 
acteristic matrix, Ea/c(< A^max), and in the data matrix 
(in the new coordinate system), Tgf^{< /Cmax), ia order 
to construct the new datacube i'^xi— ^max)- la general, 
we can discard a number of high order eigenvectors that, 
gathered together, contribute to a small fraction of the 
variance in the dataset. The method that we will use in 
order to find the number of eigenvectors we must consider 
in the analysis, which means all of the first eigenvectors 
with k < fema x: is the s cree test (see the Section see 


also Steiner et al.||2QQ9 ). It is among those eigenvectors 
limited by the scree test that we will search the presence 
of signatures of precession and rotation. Once found, 
these properties can be emphasized in the datacube re¬ 
construction process: instead of limiting the relevance as 
suggested by equation that will only discard the less 
relevant eigenvectors in the reconstruction procedure in 
order to produce a “cleaned” datacube, we can pick only 
those that contain a given physical property and pro¬ 
ceed with the r econstruction process (see the Appendix 
[b| Section Bj for more details). We will explore this in 
me Section where we will discuss reconstructed dat¬ 
acube putting in evidence properties like the jet preces- 
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Fig. 1. — Sketch of the simulation setup. In a) we show the 
general case in which we have a precessing jet observed at a given 
inclination angle 0, which is the angle between the precessing axis 
with respect to the “plane of the sky” (the x — y plane), whereas 0 
is the precessing angle of the jet axis measured from the precessing 
axis. The sense of the precession is also shown. In b) we show the 
sense of the jet rotation as seen by the observer for a rotating jet 
model in the plane of the sky. 

sion and/or rotation. 

As a final comment, the mean spectrum that has been 
subtracted from the datacube (Qy) can be reincorpo¬ 
rated in the reconstructed datacube to recover the flux- 
calibrated data. The final product is a sequence of images 
(as a function of wavelength) in the real space: 


^max))o — ^max) + Qx (9) 


from which one can reconstruct 2D images integrated in 
wavelength or select spectra from a given position in the 
images. 

We will apply this formalism to our numerical simu¬ 
lations in the next section. We use PCA to mine data 
from numerical simulations as has been done for instance 


by Heyer fc Schloerb ( 19971), [Brunt, He yer fc Mac Low| 
(2UU9|) and Carrol, F'rank & Blackman pUlUK We will 
show that rotation and/or precession can appear in inde¬ 
pendent eigenvectors, and the PCA technique can there¬ 
fore be used in the study of jet rotation to disentangle, 
in real data, signatures of rotation from other physical 
mechanisms, such as shocks and velocity asymmetries. 


3. THE NUMERICAL SIMULATIONS 

We apply the PCA analysis to the study of rotation in 
stellar jets by using synthetic emission maps generated by 
numerical simulations. The sim ulations presented her e 
are the same shown previously in Cerqueira et al. (2006). 
The simulations were perfor med with the 3D Hydro dy- 
namic code YGUAZU-a (e.g. Raga et al.||2000). Briefly, 
we run 4 models of jets with pulsation period of = 8 
years (see Table and Figure [^. The Ml model, our 
reference jet model, is intermittent, with a jet velocity 
profile given by: 




1 + Asin 



( 10 ) 


[OIF6300 [SIIF6716 


1.0 

_ 0.8 
^ 0.6 

< 0.4 
0.2 
0.0 

1.0 

_ 0.8 
^ 0.6 

< 0.4 
0.2 
0.0 

1.0 

_ 0.8 
^ 0.6 

< 0.4 
0.2 
0.0 

1.0 

_ 0.8 
^ 0.6 

< 0.4 
0.2 
0.0 

2 4 6 8 10 12 14 2 4 6 8 10 12 14 

k k 

Fig. 2.— Eigenvalues (in % of the variance) as a function of 
eigenvector’s order for models MI to M4 (from top to bottom, re¬ 
spectively), for [O l]A6300 (left) and [S Il]A67I6 (right) emission 
lines. In each panel we show also the different inclinations consid¬ 
ered: (f) = 0° (black dots), b = 15° (red dots) and (j) = 45° (blue 
dots). The ordinate has been limited to 1.2% and the first eigen¬ 
values are out of scale. However, they are given in Table where 
we also provide all values until k = 10. In each panel we have a 
horizontal, dashed black line, that was used to indicate the fcmax 
for each curve (see the text for discussion). 

to reproduce the physical conditions of the DC Tan jet. 
Besides pulsation ( model Ml), we include precession 
in M2 model, rotation in M3 model and precession and 
rotation in M4 model. For models with precession we 
adopt a period of precession equal to = 8 yr and an 
angle of precession 0 = 5°. The models with rotation 
have a rotational velocity profile given by: 

i;^ = 8kms“^^, (11) 

for R > 0.15Rj, where Rj is the jet radius {IRj = 37.4 
AU) and R is the cylindrical radius. In order to avoid 
a singularity at the jet axis, the rotational velocity is 
constant for R < O.lbRj, and equal to = 55 km s“^. 

The computational domain is a three-dimensional 
Cartesian box with spatial dimensions of x = ^ = 128 
and ^ = 512 pixels. Each pixel has a physical dimen¬ 
sion of 5.74 AU. We will present maps only for i < 128 
and j < 12^ since we are interested in the region near 
the jet inlet Che jet injection point at the computational 
domain). 

Table summarizes the main parameters of the mod¬ 
els. To l^ild the intensity matrix (equation]^ we con- 
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where vo = 300 km s ^ and A = 0.33. As discussed in 


Cerqueira et al. (2006), these values were chosen in order 


® i and j are indexes for the x and y computational cells, respec¬ 
tively. 
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-1000 100 -1000 100 -1000 100 -1000 100 

d (AU) d (AU) d (AU) d (AU) 




k = 1,[S ll]W16, (|)= 0° 
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Fig. 3.— First tomogram (T^ ^ fc, for k = 1) and eigenspec- 
trum {Ek, for k = 1) for models Ml to M4 (from left to right, 
respectively), for [O l]A6300 (two topmost panel’s sequence) and 
[S Il]A6716 (two bottommost panels) emission lines at 0 = 0°. In 
each tomogram the distance from the jet inlet, in AU, is indicated. 
The color bar indicate the level of the (normalized) intensity. The 
eigenspectrum of each model is plotted below the respective tomo¬ 
gram. 


struct velocity channel maps (VCM) around a given 
wavelength (see Cerqueira et al. 2006 for the details). 
For the forbidden emission lines |(J ljAA6300,6363 and 
[S Il]AA6716,6731p^we calculate VCM from Vrad = —400 
km s“^ to +100 km s“^, with a sampling of 10 km s“^. 
Each channel map is one of the m slices in A of the dat- 
acube, and we can build the intensity matrix {Iijx)o foi* 
each A (or, alternatively, m or Vrad)- 
In order to consider the effect of finite angular resolu¬ 
tion we convolved each VCM with a Gaussian PSF with 
FWHM of ~ 5 pixels (^ 28.7 A U), or 0721 for a jet at 


a distance of Ta urus (140 pc; see Kenyon, Dobrzycka & 
Hartmann||1994 ), for example. 


4. EIGENVALUES, TOMOGRAMS AND EIGENSPEGTRA 
4.1. Eigenvalues and the seree test. 

In the PC A, it is important to know how many eigen¬ 
vectors one must consider in the analysis. One way to 
deter mine this is to apply the seree test (Steiner et al. 
2009), which is a visual test that is used to find the order 
k where the eigenvalue levels off. To look at these /c’s 
for our different models, emission lines and inclination 
angles with respect to the plane of the sky, we built di¬ 
agrams of eigenvalues A/, as a function of k. In Figure 


As we have already mentioned and discussed, results will be 
presented only for the [O l]A6300 and for the [S Il]A6716 emission 
lines. 


Fig. 4.— The same as in Figure but for k = 2 (the second 
tomogram, for = 2, and its respective eigenspectrum, 

for k = 2). 

[^we show the eigenvalues for models Ml, M2, M3 and 
M4 (from top to bottom, respectively), for [O IJA6300 
(left) and [S IIJA6716 (right). The dots of different col¬ 
ors are associated with a specific inclination angle: black 
for 0 = 0°; red for 0 = 15°; blue for 0 = 45°. We 
note from Figure that the eigenvalues go to ^ 0 after a 
given k = /Cmax- In order to help to find these /Cmax, we 
have drawn a horizontal line in each plot at a constant 
Athreshoid = 10“^ %. The first point (for each model, in¬ 
clination and emission line) intercepted by this horizontal 
line defines the /Cmax- In Appendix |A] we present in Ta¬ 
ble [| the first ten eigenvalues for eadi one of the models. 
We present also the number of relevant eigenvectors as 
determined by the criteria of the scree test (see Table 
1^. This means that for each combination of parameters 
(model, emission line and inclination angle) we have an 
already predefined number of eigenvectors that we must 
look in searching for relevant uncorrelated phenomena. 
However, we used higher order eigenvectors to perform a 
complete analysis of the problem, although they will not 
be presented here due to the lack of space. 


4.2. Tomograms and eigenspeetra 

To obtain an image in the new coo rdinate system, or 
tomo gram (or also an eigenimage; see |Heyer fc Schloerb 


1997), we use equations 
data matrix (see equationi 


and A column in the 
can be then transformed in 


a image, that is the projection of the datacube (in the 
new coordinate system) onto the chosen eigenvector (of 
order k). The eigenvector, on the other hand, is a linear 
combination of the original (spectral) coordinates, and 
its coefficients, or weights, can be both positive or nega- 
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k = 2, [O l]X6300, (^=15 
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k = 2, [S ll]W16, (|)=15° 




Fig. 5.— Second tomogram (T^ ^ for k = 2) and eigenspec- 
trum {Ek, for k = 2) for models Ml to M4 (from left to right, 
respectively), for [O l]A6300 (two topmost panel’s sequence) and 
[S Il]A6716 (two bottommost panels) emission lines at (/) = 15°. In 
each tomogram the distance from the jet inlet, in AU, is indicated. 
The color bar indicate the level of the (normalized) intensity. The 
eigenspectrum of each model is plotted below the respective tomo¬ 
gram. 

tivj^ Viewed as a plot of weight versus radial velocity, 
the eigenvector is called eigenspectrum. Therefore, a to¬ 
mogram can displays either positive or negative regions^ 
that are related with the weights in the corresponding 
eigenspectrum. This will allow us to investigate both 
and identify some kinematical “properties”, like jet pre¬ 
cession or rotation. 

4.2.1. Jet or precessing axes in the plane of the sky (cj) = 0°) 

In Figure]^ we show the first tomogram {k = 1) and its 
associated eigenspectrum for models Ml to M4 (from left 
to right), for [O I]A6300 (top) and [S II]A6716 (bottom) 
emission lines, at 0 = 0°. The color bar in the tomo¬ 
grams (from here and after) represent a non-dimensional 
quantity, since pixel’s intensities have been normalized 
to the maximum value in each tomogram. Distances, in 
AU, are indicated in the ordinate of the leftmost tomo¬ 
gram, as well as in each abscissa. For the eigenspectra, 
the weights are plotted as a function of the radial veloc¬ 
ity in km s“^. These k = 1 eigenvectors for each model 
(see Table present tomograms that are similar to an 
integrated image in the respective emission line. They 
emphasize the presence of the two working surfaces: one 
internal near to the jet inlet, and the leading one at the 
tip of the jet, where most of the emission comes from 

The only constraint is that their sum in quadrature must be 
equal to the unity, because the eigenvectors are normalized. 



Fig. 6. — a) Geometry for a jet inclined towards the observer. 
Lateral backflowing material with velocities compose with 

the jet velocity and favors high radial velocities in the half bottom 
of the jet cross section (as seen by the observer; indicated by the 
blue arrow and with velocity Vr2) in comparison with the half top 
part of the jet cross section (red arrow with velocity Vri)- Both are 
blueshifted velocities but appear in the tomogram of the model Ml 
in of Figure [2 as a red and blue regions, b) An inclined bow shock 
can also expKin the tomograms/eigenspectra of the model Ml and 
cannot be ruled out in the case of the internal working surface, 
where the internal bow shock propagates faster than the leading 
one due to the non-stationarity of the ambient medium ahead of 
it. 

the radiative losses behind the shocks. This eigenvector 
traces the most important shocks in the system. The 
respective eigenspectra peak at ^ 0 km s“^ in the case 
of models Ml and M3, and at ^ 10 km s“^ in the pre¬ 
cessing cases M2 and M4 (M4 also has rotation besides 
precession). 

In Figure we show the second tomogram {k = 2) 
and its associated eigenspectrum for models Ml to M4 
(from left to right), for [0 l]A6300 (top) and [S II]A6716 
(bottom) emission lines, at 0 = 0°. 

For Ml model (panels in the first column), this eigen¬ 
vector contributes to 0.2597% and 0.2655% of the vari¬ 
ance (for [O I] and [S II], respectively; see Table and 
traces velocity gradients at the jet head (for bothemis- 
sion lines; their tomograms and eigenspectra are equiva¬ 
lent) that are perpendicular to the jet axis, but symmet¬ 
ric with respect to it. 

The tomograms for the M2 (precessing) model (panels 
in the second column of the Figure ffl are considerably 
different for the two emission lines. In the case of the 
[O I]A6300, the tomogram combined with its eigenspec¬ 
trum clearly traces the precession of the jet axis. The 
eigenspectrum has anti correlated wings: the red wing 
of the line, associated with positive weights and posi¬ 
tive regions in the tomogram, traces the receding portion 
of the jet (coincident with the internal working surface 
near to the jet inlet), while the blue wing of the line, 
associated with negative weights and negative regions in 
the tomogram, traces the approaching jet. We note that 
there is also a positive gradient, from the internal work¬ 
ing surface (hereafter, IWS) up to the jet head, in the 
tomogram’s intensity along the jet. The eigenspectrum 
for the [S II]A6716 emission line also depicts the same 
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Fig. 7.— The same as in Figure but for k = 3 (the third 
tomogram, for k = 3, and its respective eigenspectrum, 

for k = 3). 

anti correlation in the wings of the line, but here the in¬ 
terpretation is not as straightforward as before. In the 
eigenspectrum, there is only one peak for each positive 
and negative weights. This suggests that all blue regions 
in the tomograms are blueshifted with respect to the red 
ones. According to this interpretation, the IWS is reced¬ 
ing and there is a smooth gradient in blue regions from 
the IWS towards the jet head suggesting the precession. 
However, at the jet head seen in the tomogram we have 
a strong gradient in radial velocity, and the anti corre¬ 
lation seen in the eigenspectra in this case may be due 
to a combination of these two features, and could not be 
attributed only to the precession. We note that a sig¬ 
nature for the precession can be found in all tomograms 
until k = 4 (not shown here). 

The second eigenvector for M3 (rotating) model ac¬ 
counts for 6.7645% and 3.9785% of the variance in the 
dataset for [O I]A6300 and [S II]A6716 emission lines, 
respectively (Table [^. Their tomogram/eigenspectrum 
are in the third column of the Figure]^ The eigenspec¬ 
trum for the [O I]A6300 emission line indicates again an 
anti correlation between both wings of the line. Nega¬ 
tive weights are associated with negative regions (blue 
regions) in the tomogram, indicating a receding jet (it 
is in the red wing of the line), while positive weights, in 
the blue wing of the line, indicate that the red regions in 
the tomogram are approaching. This is consistent with 
the sense of the rotation imposed at the beginning of our 
simulation. This eigenvector for [O I] traces, then, the 
rotation of the jet. Conversely, the same tomogram and 
the eigenspectrum for the [S II]A6716 emission line are 
completely symmetric with respect to the jet axis and to 


Fig. 8. — Third toruograui (T^ ^ fc, for k = 3) aud eigeuspec- 
truru {Ek, for k = 3) for ruodels Ml to M4 (froui left to right, 
respectively), for [O l]A6300 (two topuiost pauel’s sequeuce) aud 
[S Il]A6716 (two bottouiuiost pauels) euiissiou hues at 0 = 45°. lu 
each toruograui the distauce froru the jet iulet, iu AU, is iudicated. 
The color bar iudicate the level of the (uorrualized) iuteusity. The 
eigeuspectruru of each ruodel is plotted below the respective toruo- 
graru. 

the zero radial velocity, respectively. There is no trace 
at all for the presence of the rotation in this tomogram 
for this emission line. It will appear only in the third 
eigenvector (not show here). 

The second eigenvector for M4 (rotating and pro¬ 
cessing) model accounts for 18.007% and 12.7514% 
of the variance of the data, for the [O I]A6300 and 
[S II]A6716 emission lines, respectively. Their tomo- 
gramVeigenspectrum are in the fourth column of the Fig¬ 
ure The tomogram and the eigenspectrum for the 
[O IJA6300 emission line shows anti correlation between 
both sides of the jet axis near to the jet inlet only. Posi¬ 
tive regions, red in the tomogram, are correlated with the 
positive, red wing of the line, while negative regions, blue 
in the tomogram, are correlated with negative weights in 
the blue wing of the line. This is consistent with the 
interpretation that the jet is rotating, and the correct 
jet rotation sense is recovered. From the IWS until the 
jet head the jet is dominated by blue, negative regions. 
The gradient seen in the tomogram can suggest a pre¬ 
cession. For the [S II]A6716 emission line the situation is 
less clear. Some asymmetries with respect to the jet axis 
can be detected near to the jet inlet, but it is not obvi¬ 
ous that rotation is involved. The eigenspectrum seems 
to reflects the gradient in radial velocity at the jet head, 
which is in turn very similar to the one detected in the 
M2 model for this same tomogram. 

In the following sections, we will not show the first 
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Fig. 9.— The same as in Figure but for k = 4 (the fourth 
tomogram, for k = 4, and its respective eigenspectrum, 

for k = 4). 

tomograms anymore (i.e., tomograms for k = 1): they 
are always similar to an image in the real space of the 
integrated emission line and they do not contribute to 
our present purposes, which is the search for evidence of 
rotation and precession in jets. 


4.2.2. Tomograms for moderately inelined systems (cf = 15°^ 

In Figurewe show the second tomogram (k = 2) and 
its associatea eigenspectrum for models Ml to M4 (from 
left to right), for [O I]A6300 (top) and [S II]A6716 (bot¬ 
tom) emission lines, at 0 = 15°. All eigenspectra are 
blueshifted in comparison with the models in previous 
section, as expected for jets that are pointing towards the 
observer. They also present anti correlation between the 
red and blue wings of the line, no matter if we are look¬ 
ing for a simple, non-precessing and non-rotating (Ml) 
model or a complex (M4) one, and this fact deserves a 
careful analysis. 

In the case of the reference model Ml, the structures 
(positive/negative regions) in the tomogram are symmet¬ 
ric with respect to the jet axis. The portions of the jet 
immediately behind the working surfaces (at d{y) <150 
AU and at 550 AU < d{y) < 600 AU) are correlated with 
the negative weights of the eigenspectra for both emis¬ 
sion lines (see the first column of Figure |^, and then, 
they are correlated with the blue wing ofthe emission 
line. Beyond these limits (that is, for d{y) ^ 180 AU 
and d{y) ^ 630 AU), the tomograms display positives 
(see the color bar in each figure) regions, that are corre¬ 
lated with the red wing of the line. We interpret these 
structures at the jet working surfaces as a signature of the 
backflowing post-shocked material, that escapes laterally 






Fig. 10.— luiages (left) of the collapsed origiual datacube for the 
M3 (rotatiug) uiodel arouud the [O l]A6300 (top) aud [S Il]A6716 
euiissiou hues. Spectra (right) have beeu extracted froru the area 
iudicated by the white circles iu the iuiages usiug the origiual dat¬ 
acube (black solid hue), aud the recoustructed datacubes cousid- 
eriug eigeuvectors uutil fcmax = 1 (red), Aimax = 2 (blue), fcmax = 3 
(greeu) aud fcmax = 4 (black stars). Fluxes are uorrualized to uuity, 
aud each ruaxiuiuru is iudicated iu the right side of the color bars. 
Iu the extracted spectra, iuteusities are also uorrualized with re¬ 
spect to the ruaxiuiuru of the iutegrated spectra (black solid hue). 
The circles iu the iuiages have a diauieter of 4 pixels, or 22.9 AU. 
This correspouds to 07l6 cousideriug the distauce to Taurus (140 
pc), for exauiple. The origiu of the coordiuate systeiu correspouds 
to the positiou of the jet iulet, aud distauces froru the origiu are 
iudicated iu AU. 


at the Mack disk (or jet shock), favoring lower (absolute) 
radial velocities on the top half of the jet cross section, 
and, analogously, higher (absolute) radial velocities on 
the bottom half of the jet cross section. In Figure (a) 
we present a cartoon summarizing the interpretation. We 
note that an inclined bow shock could also account for 
the blue/red asymmetry in the line and cannot be ruled 
out (see the sketch presented also in Figure [^b), partic¬ 
ularly for the case of an IWS, for which the bow shock 
propagates into a non-stationary ambient medium. In 
this case this effect can be more noticed since we expect 
high bow shock velocities in comparison with the leading 
bow shock. 

The precessing model M2 also show an asymmetric 
eigenspectrum. In the [O IJA6300 emission line, the to¬ 
mogram essentially traces the IWS near the jet inlet. It 
is a positive region (red in the color bar), correlated with 
the positive weights that are present in the red wing of 
the line, as we can see in the eigenspectrum. There is a 
clear contrast between positive (internal emission knot) 
and negative regions (the jet head), suggesting a preces¬ 
sion. However, positive weights dominate the eigenspec¬ 
trum (it has a maximum that is more than three times 
the negative peak at Vrad = ~90 km s“^), which is evi¬ 
dent also in the tomogram. For the [S IIJA6716 emission 
line, the contrast between positive and negative weights 
is less intense and the structures more clear. There is a 
gradient in the negative region that goes from the IWS 
and peaks at the jet head. There, we have also a pos¬ 
itive region, suggesting that the same effect explored in 
Figure may also be present. However, the precession 
breaks the symmetry with respect to the jet axis. 
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Fig. 11.— Images (left) and spectra (right) from the recon¬ 
structed datacube for the precessing model M2 at 0 = 0°, for the 
[O l]A6300 (top) and [S Il]A6716 (bottom) emission lines. Eigen¬ 
vectors from 2 to 4 were used in the reconstruction process, since 
they have signatures for the jet precession. Circles labeled from 
A to E mark the positions where a spectrum has been extracted. 
They are in the right side panels, normalized to unit. We added a 
constant 2n, where n = 0 —4 to the spectra at positions A to E, re¬ 
spectively, in order to separate them in different levels, to improve 
the visualization. The images are integrated in wavelength. 


For the M3 (rotating) model (third column in Figure 
1^, the second tomogram shows positive (red) and nega¬ 
tive (blue) regions that are distorted with respect to the 
jet axis. The region that extends from the jet inlet until 
the IWS is blue in the left side of the jet axis for both 
emission lines (although more clearly seen in the tomo¬ 
gram associated with the [O IJA6300 emission line). At 
the tip of the IWS, the red (positive) region is distorted 
toward the right side of the jet axis, indicating that this 
side of the jet is more redshifted than the other one (since 
positive weights in the eigenspectra are in the red wing 
of the line; see the eigenspectrum associated with this 
model in the third column of Figure]^. At the jet head, 
the rotation signature is evident only in the case of the 
[S IIJA6716 emission line: we can interpret the blue/red 
region in the jet head in this tomogram as in the case of 
model Ml, and attributes to the rotation the cause for 
the observed asymmetry with respect to the jet axis. 

The second eigenvector for the model M4 contributes 
to 33.2961% and 14.5136% of the variance in the dataset 
for [O IJA6300 and [S IIJA6716 emission lines, respec¬ 
tively, at 0 = 15°. The tomogram and eigenspectrum 
can be seen in the last column of the Figure Again, 
the eigenspectra for both emission lines show blue/red 
wing asymmetry. The tomograms are, however, differ¬ 
ent. For the [O IJA6300 emission line, positive regions in 
the tomogram are concentrated in the IWS. It has a peak 
on the right side of the jet axis, showing an asymmetric 
distribution in intensity with respect to the jet axis (as 
in the case of the M3 model). All other parts of the jet 
are mostly dominated by blue, negative regions, with a 
maximum at the jet head. There is a negative gradient 
from the jet inlet to the IWS, and then a positive gradi¬ 
ent from the IWS up to the jet head, which is suggest the 
precession. Although the eigenspectra are similar for the 


[S IIJA6716 emission line in comparison with the one for 
the [O IJA6300 emission line, their tomograms are not. 
Both positive and negative regions peaks at the jet head, 
suggesting that the anti correlation in the wings of the 
line are related to gradients in velocity at the jet region. 

In Figure we show the tomograms (first and third 
rows of panels) and eigenspectra (second and fourth rows 
of panels) for models Ml to M4 (from left to right) at 

= 15° and k = 3. As before, the two topmost pan¬ 
els refer to the [O IJA6300 emission line, while the two 
bottommost panels concern the [S IIJA6716 emission line. 
There is a pattern in the eigenspectra that is almost the 
same for all case a^ there is an interval in radial velocity 
of negative weights, bracketed by two peaks of positive 
weights, which make the interpretation far less obvious. 
In particular, a negative region in a given tomogram can 
be either red- or blue-shifted with respect to a positive 
one. The same “color” in a tomogram can, in this case, 
be indicative of different kinematics, and the interpreta¬ 
tion is actually difficult. 

Ml model (leftmost column in Figure has tomo¬ 
grams that show negative/positive regions that peak at 
the jet head (for both emission lines). The positive (red) 
regions at the jet head are correlated now with positive 
weights in the eigenspectra, since the maximum in both 
tomograms is localized in this region. Then, the negative 
region just above it in the case of [O IJA6300 emission 
line, and the negative regions symmetrically displaced 
with respect to the jet axis in the case of [S IIJA6716 
emission line are redshifted in comparison with the red 
regions at the jet head, as we can see in the eigenspec¬ 
trum. This interpretation is still compatible with the 
scenario proposed in Figure The same is true for the 
IWS, which correlates with the positive/negative peaks 
at Vrad = “50 and -80 km s“^, respectively. The yellow 
(positive) region in the jet base can be only correlated 
with the positive weights at the bluemost wing of the line 
in the eigenspectra. 

M2 model (panels in the second column of Figure 
shows that the highest values for the positive and neg¬ 
ative regions in the tomogram occur at the jet head. 
Negative regions are concentrated at the jet tip, while 
the positive regions are immediately behind of it, with 
a negative gradient of intensities from the jet head until 
the IWS. The peak in the positive weights of both eigen¬ 
spectra (i.e., for both emission lines) is at Vrad = —100 
km s“^, that are in the blue wing of the lines. This 
means that the negative blue region on top of these pos¬ 
itive regions at the jet head are redshifted with respect 
to them, since the negative weights in the eigenspectra 
are at lower (absolute) radial velocities (^ —60 km s“^). 

The third tomogram for the M3 and M4 models (third 
and fourth columns of panels in Figure respectively) 
at (j) = 15° show a twisted pattern in the regions of 
positive and negative values, which is more evident for 
the [O IJA6300 emission line in comparison with the [S 
IIJA6716 emission lin^^ We can interpret this behav¬ 
ior saying that negative regions in the tomograms can 

The only exception is the model M4, in the [S IIJA6716 emis¬ 
sion line (see below) 

It becomes also evident in high order tomograms for the [S 
II] A6716 emission line. We should note that it it also present in high 
order tomograms for the non-inclined, 0 = 0°, rotating models. 
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Fig. 12.— Left:) The inclination of the jet axis with respect to the plane of the sky, here defined as z, as a function of the vertical distance 
from the jet inlet for the model M2 (precessing jet), estimated from the simulation (solid line), using velocities taken in a straight line along 
the y—coordinate, starting at the jet inlet at x = y = 0. The dashed line is an arbitrary cosine function. Middle:) Normalized fiuxes taken 
from the reconstructed datacube considering only the second eigenvector (k = 2) (solid line: [O l]A6300; dashed line: [S Il]A6716). Right:) 
The same as in the middle panel but using the eigenvectors 2 < A: < 4 to reconstruct the datacube. The fiuxes have been normalized to 
their maximum value in each curve. In the rightmost panel, we have limited the ordinate to 0.2. 


be both redshifted or blueshifted with respect to posi¬ 
tive regions placed in a symmetric position with respect 
to the jet axis. This is because the negative peaks in 
the eigenspectra for both models are bracketed by two 
positive peaks, and that it is an effect of the rotation, 
that changes the tomograms with respect to the non¬ 
rotating Ml and M2 models even if the eigenspectra have 
the same shape/form. However, by analyzing only this 
eigenvector, the presence of such a twisted feature in the 
tomograms does not allow us to say/conclude anything 
about the rotation sense of the jet. 

Unlike the case of 0 = 0° shown in the previous sec¬ 
tion, evidences for the rotation here are more subtle. The 
mere fact that we are “observing” the jet at a non-null 
inclination angle makes the analysis of the principal com¬ 
ponents more difficult. 

4.2.3. Tomograms for highly inclined systems ((f — 45°/ 

At 6 = 45°, the dataset becomes even less redun- 
danip^ with an appreciable amount of variance dis¬ 
tributed in higher order eigenvectors. In Figures and 

we show as an example the third and fourth tomo- 
grams / eigenspectrqr respectively, for models Ml to M4 
at 0 = 45°. As b^re, the results for the [O IJA6300 
(topmost two panels) and [S IIJA6716 (bottommost two 
panels) emission lines are shown. Some of the patterns 
already discussed can be seen here too. Eigenspectra 
becomes more complex, in part due to gradients in ve¬ 
locity at the jet head. The clumpy, alternating pattern 
observed there specially in models M3 and M4 (for k = 3 
and k = 4) makes the kinematic analysis, and its inter¬ 
pretation, more difficult. The presence of the rotation, 
for instance, is only suggested near the jet inlet in the 
k = A tomograms of model M3 (and it is not so evident 
for the M4 model; see Figure 1^. It is present also in 
higher order tomograms (not shown here). 

In the sense that the number of relevant eigenvectors, defined 
by the scree test for example, is higher in comparison with a more 
redundant dataset. In other words: if just a few eigenvectors can 
be used to describe the dataset and the majority of them can be 
discarded for this purpose, then it is a highly redundant dataset. 

The first two tomograms in these cases traces mainly the IWS 
and the leading working surfaces and will not be discussed here. 


We conclude that it is not possible to express, or to de¬ 
scribe, the precession and/or the rotation with a single 
eigenvector, unless if the precession/jet axis lies in the 
plane of the sky. We have actually a collection of eigen¬ 
vectors for which the signature for the precession and/or 
for the rotation can be found. In general, the higher the 
inclination angle 0, the higher will be the eigenvector’s 
order in which they will manifest. In this Section we have 
discussed the most relevant eigenvectors. Furthermore, 
the interpretation of higher order tomograms becomes 
less and less evident, with eigenspectra showing several 
negative and positive successive peaks. 

5. PCA SIGNATURES FOR PRECESSION AND ROTATION 
FROM RECONSTRUCTED DATACUBES 

In this section we will present images from recon¬ 
structed datacubes, obtained using the equation ^ for 
a collection of eigenvectors that are representative of a 
given feature (see Appendix [B|). We will also discuss the 
properties of spectra extractM at fixed spatial positions 
obtained from reconstructed datacubes. We note that all 
images presented in this section are actually integrated 
in wavelength. 

Before entering into the details of our results, and as 
a test for the method, we can choose a region in a given 
image and reconstruct, step by step, the line profile us¬ 
ing a certain number of eigenvectors. In Figure we 
show the integrated (non treated by the PCA) images 
of the [O I]A6300 (top panel) and [S II]A67I6 emission 
lines (bottom panel) for the M3 (rotating) model. Su¬ 
perimposed on the images we show circles (of 22.9 AU of 
diameters; which corresponds to 07l6 at 140 pc of dis¬ 
tance) to indicate the position where spectra have been 
extracted. The spectra are on the right panels (black, 
solid lines). We also show line profiles obtained from 
the reconstructed datacube, keeping in the reconstruc¬ 
tion procedure an increasing number of eigenvectors from 
I to 4 (see equation |^: /cmax = 1 (red curve), /cmax = 2 
(blue curve), /Cmax (green curve) and /Cmax = 4 (black 
stars). In this example, the reconstructed profiles con¬ 
verge very well to the original one for both emission lines 
(the properties and features of the profiles will be dis¬ 
cussed in the next Section), using only 4 eigenvectors. 
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Fig. 13.— Reconstructed images (left) and extracted spectra 
(right) for Model M3 at 0 = 0° and for [O l]A6300 (top) and 
[S Il]A6716 (bottom) emission lines. Different eigenvectors have 
been considered to reconstruct the datacube: k = 2 and 4 for 
[O l]A6300 and k = 3 and 4 for [S Il]A6716 lines, as indicated 
in the bottom left part of each image. Spectra were extracted at 
symmetrical positions, labeled A, B and C, with respect to the jet 
axis. The colors of the solid lines (right panels) corresponds to 
the color of the slits. Distances from the jet inlet are indicated in 
AU in images. Spectra are normalized and we have added to their 
intensity a constant 3n, with n = 0, 1 and 2 at the positions A, B 
and C, respectively. The images are integrated in wavelength. 


5.1. Jets in the plane of the sky 

In Figure pT] we show the image (left panels) and spec¬ 
tra (right panels) obtained from the reconstructed dat¬ 
acube for the model M2 (precessing model) at 0 = 0° 
(which means that the axis of precession is coincident 
with the axis). As discussed in Section (Section 
|4.2.1 ), we have identified the precession in eigenvectors’ 
orders 2 < k < 4 amon^he eigenvectors suggested by 
the scree test (see Table^. The datacube was then re¬ 
constructed considering these three eigenvector^^ For 
the sake of clarity and to allow the comparison between 
them, the results for the [O I]A6300 and [S II]A6716 emis¬ 
sion lines are shown. The emission line is indicated in the 
top-left corner of each image. We choose different posi¬ 
tions along the jet to extract a spectrum. These regions 
are indicated by circles superimposed in the image (as 
in Figure the size of the circle is representative of 
the collected area which means 07l6 at a distance of 140 
pc, for example). Viewing from the base of the jet (re¬ 
gion A) up to its tip (region E), each spectrum peaks at 
different radial velocities, suggesting clearly a precession 
pattern. Assuming as an upper limit for the jet velocity 
< 400 km s“^ (see equation 10), we can estimate 
the radial velocity to be of the order oi|'r’^ad| < '^i,zimsin^, 
or I'T’radl ^ 34 km s“^. The emission peaks fall well 
within this limit. 

Figure indicates that there is a variation in the 
flux along the jet length that can be due to changes 
of the viewing angle. In Figure 12 we try to correlate 


At each image obtained from a reconstruction process we in¬ 
dicate inside the parenthesis of the expression placed at the 

left-bottom corner of the images, the selected eigenvectors used in 
the reconstruction process. 


Fig. 14.— Reconstructed images (left) and extracted spectra 
(right) for Model M4 at 0 = 0° and for [O l]A6300 (top) and 
[S II]A6716 (bottom) emission lines. Different eigenvectors have 
been considered to reconstruct the datacube: = 2 to 8 for the 

oxygen and A: = 2 to 7 for the sulfur lines, as indicated in the bot¬ 
tom left part of each image. Spectra were extracted at symmetrical 
positions, labeled A, B and C, with respect to the jet axis. The 
colors of the solid lines (right panels) corresponds to the color of 
the slits. Distances from the jet inlet are indicated in AU. Spectra 
are normalized and we have added to their intensity a constant 
3n, with n = 0, 1 and 2 for spectra at the positions A, B and C, 
respectively. The images are integrated in wavelength. 


these quantities. For a fixed x—coordinate at the jet in¬ 
let {x = 0), we compute tan“^(t;^ad/'^?/) along the jet 
length, which is an estimate for the angle i that the jet 
axis makes with the precession axis that lies in the plane 
of the sky (the y axis at x = 0). In this equation, Vrad 
and Vy are the velocities parallel and perpendicular to 
the line of sight, respectively, taken from the simulation. 
This quantity is plotted in the leftmost panel in Figure 


12 (solid line), where we show also an arbitrary cosine 
Tunction (dashed line). In the middle and right pan¬ 
els we show the flux in the [O I]A6300 (solid line) and 
[S II]A6716 (dashed line) emission lines, taken along the 
precession axis from the reconstructed datacube, consid¬ 
ering the eigenvectors with k = 2 (middle panel) and 
2 < /c < 4, respectively. The abrupt variation of the 
intensity along the jet seems to correlates with the vari¬ 
ation in the inclination angle when we consider the com¬ 
bined eigenvectors (rightmost panelj^ 

In Figure we show images (left) and spectra (right) 
from the reconstructed datacube for the [O I]A6300 (top) 
and [S II]A6716 (bottom) emission lines, for the rotating 
M3 model observed at 0 = 0°. For the [O I]A6300 line, 
we have found evidence for rotation in the eigenvectors 
2, 4, 7, 9, 12 and 14 (this is not a comprehensive list). 
For the [S II]A6716 line, the eigenvectors associated with 
the rotation are 3, 4, 7, 9, 12 and 14. The scree test (see 
Section suggests for this model, however, that the rel¬ 
evant e^envectors are those with k < Q (see Table 1^ and 
Figure^. We have, then, restricted the relevance m the 
eigenvectors in the reconstruction procedure up to this 


The mean spectra, Qx, has not been summed up in the re¬ 
construction process, since we took just a few eigenvectors to re¬ 
construct the datacube. 
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Fig. 15.— Reconstructed images (left) and spectra (right) for 
the precessing model M2 for 0 = 15° (first and second panels, 
from top to bottom) and 0 = 45° (last two bottom panels), for [O 
1JA6300 (first and third rows) and [S 11JA6716 (second and fourth 
rows) emission lines. Different eigenvectors have been used in the 
reconstruction procedure. Their range is show at the top left of 
each image. Circles labeled from A to E mark the positions where 
a spectrum has been extracted. They are in the right side panel, 
normalized to unit (we added a constant 2n, n = 0 — 4 from A 
to E, in each spectra in order to separate them in different levels). 
The images are integrated in wavelength. 


limit. In these two images we placed artificial slits at 
increasing distances from the jet inlet, at positions la¬ 
beled A, B and C. The slits are symmetrically displaced 
in pairs with respect to the jet axis. The extracted spec¬ 
tra are plotted in the right panels (the colors of solid lines 
used to draw each spectra are associated with the color 
of the “apertures” in the FoV). The peaks of red/black 
curves suggest that the jet material at the right side of 
the jet axis is receding, while the left side is approaching, 
recovering the sense of the jet rotation imposed initially, 
independent of the position along the jet axis A, B or C. 

In Figure 14 we show images (left) and spectra 
(right) from tne reconstructed datacube for the (rotat- 
ing/precessing) model M4, in the [O I]A6300 (top) and 
[S II]A67I6 emission lines. For this case, we have evi- 
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Fig. 16.— Reconstructed images (left) and extracted spectra 
(right) for Model M3 at 0 = 15° and for [O l]A6300 (top) and 
[S II]A6716 (bottom) emission lines. To reconstruct the datacube 
we have used the eigenvectors A: = 2, 3, 4 and 6 for oxygen and 
A; = 2 to 5 for the sulfur emission lines, as indicated in the bottom 
left part of each image. Spectra were extracted at symmetrical 
positions with respect to the jet axis, labeled A, B and C. The 
colors of the solid lines (right panels) corresponds to the color of 
the slits. Distances from the jet inlet are indicated in AU. Spectra 
are normalized and we have added to their intensity a constant 3n, 
with n = 0, 1 and 2 at the positions A, B and C, respectively. The 
images are integrated in wavelength. 

deuce for rotation and precession in several eigenvectors 
with order higher than 2. We then keep the eigenvectors 
from 2 to 8 and from 2 to 7 in the reconstruction process 
of the datacube for oxygen and sulfur emission lines, re¬ 
spectively, obeying the limits suggested by the scree test 
(see Table 1^. We have extracted spectra in pairs that do 
not follow a straight line along the axis but instead, 
we place the slits along a “suggested” jet symmetry axis 
(see Figure fT3. We see in the spectra that there is ev¬ 
idence for the precession: radial velocities peaks for slit 
pairs changes from positive to negative values from B to 
C, indicating that jet material in B is receding while in 
C it is approaching, as expected. There is also an indi¬ 
cation for rotation, since the red curves peaks are always 
redshifted with respect to the black curve, as we should 
expect considering the sense of the jet rotatiorj^ 


5.2. Inclined systems 

In Figure \lb\ we show the images (left) and spectra 
(right) for the (precessing) model M2 at 0 = 15° (two 
upper panels) and at 0 = 45° (two bottom panels), for 
different emission lines, as indicated in each figure. Each 
image/spectrum has been obtained from a reconstructed 
datacube, using those eigenvectors that indicate the pres¬ 
ence of the precession in each case: /c = 2, 3, 4 and 6 for 
the 0 = 15° (first and second rows); 2 < /c < 7 for 
the oxygen line at 0 = 45° (third row of panels) and 
2 < /c < 5 for sulfur line at 0 = 45° (bottom panel). 
We can see that the peaks of the emission lines indicate 


We have also placed the slits parallel to the y-axis, and take 
into account the lateral shifts in x in order to “follow” the jet 
precession. The results are qualitatively similar to that presented 
in Figure the only difference being a small shift in the radial 
velocity at position A. 
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Fig. 17.— Images (left) and spectra (right) for model M4 at 
(j) = 45°, for [O l]A6300 (top) and [S Il]A6716 (bottom) emission 
lines. In the reconstruction process of the datacube we have used 
eigenvectors = 4, 5, 7, 8, 11, 12, 13 and 14 for the oxygen emission 
line and = 2 to 6, and k = 9 and 11 for the sulfur emission 

line, as indicated in the bottom left part of each image. Spectra 
were extracted in pairs, along and across the jet length, at different 
positions labeled A, B and C. The colors of the solid lines (right 
panels) corresponds to the color of the slits. Distances from the 
jet inlet are indicated in AU in images. Spectra are normalized 
and we have added to their intensity a constant 3n, with n = 0, 1 
and 2 at the positions A, B and C, respectively. The images are 
integrated in wavelength. 


the presence of precession, and that it is consistent with 
the pattern of the precession obtained for model M2 (for 
0 = 0°; left panel in Figure 12). As we increase the 
inclination angle, the amplitude of the variation of the 
radial velocity increases (bottom panels). The scenario 
is independent of the emission line, although maps from 
[O I] and [S II] may differ considerably. 

Eigenvectors 2, 3, 4 and 6 have been used to recon¬ 
struct the datacube for the (rotating) M3 model, at 
0 = 15° for the [O IJA6300 emission line, while eigenvec¬ 
tors from 2 to 5 were used to reconstruct the datacube 
around the sulfur line for the same model and inclination 
angle. In Figure pT] we show, as in the previous figures, 
the reconstructednnage and several spectra, extracted 
along and across the jet axis. As in the case of 0 = 0° 
(see Figure spectra from red slits (right side of the 
jet axis) aremways redshifted in comparison with spec¬ 
tra from black slits (positioned on the left side of the jet 
axis), suggesting that the jet is rotating and recovering 
the initial rotating sense. These peaks occur, however, 
at different radial velocities. While in A and C the spec¬ 
tra peaks at Vrad ^ —100 km s“^, in B they peaks at 
'^rad ^ “00 km s“^. This difference may be due to the 
presence of an internal working surface in B. Another 
potential source for this difference is the fact that jet ve¬ 
locity is time variable (see equation 10), and this lower 
(absolute) value for the radial velocity is consistent with 
the one expected when the jet velocity is in its minimum, 
'^j,min = 200 km s“^, when projected at 0 = 15°. 

In the more complex case in which we have rotation 
and precession, as in model M4, we can see the same 
effect (in all the emission lines). In Figure 17 we show 
images (left) and spectra (right) for the moHS M4, for 


Fig. 18.— Array of slits superimposed on the reconstructed im¬ 
age, zoomed at the jet inlet (left panels). The array of 4x7 slits 
disposed across and al ong the jet length is equivalent to that of 
|Cerqueira et al.|(|2006|. The open circles in the right panels are 
tFe mean values oTThe radial velocity shift from symmetrically 
disposed slits (3-5), (2-6) and (1-7), considering the four differ¬ 
ent regions (from I to IV). The error bars indicate the dispersion 
around the mean value. The red solid line is obtained calculating 
analytically Av^ad = 2 x Vradi where v^ad = Scoscp • Rj/R. In 
the top panels we have the result for M3 model at 0 = 0° and in 
the bottom panel, 0 = 15°. In all cases, the emission line is the 
[O l]A6300. The images are integrated in wavelength. 


[O I]A6300 (top) and [S II]A6716 emission lines. The 
eigenvectors used in the reconstruction process are in¬ 
dicated in the bottom-left corner of each figure. The 
pattern observed in the spectra is consistent with a ro¬ 
tating and velocity variable jet. However, the spectra 
are substantially more fuzzy in comparison with those 
obtained for this same model and emission lines at small 
inclination angles (see Figure 14). 

In Figure we attempt to recover the initially im¬ 
posed rotational profile using the reconstructed dat¬ 
acube. We build an array of 4x7 slits, or four regions 
(I to IV) at increasing distances from the jet inlet and 
seven positions (1 to 7) across the jet axis. The size of the 
slit s and their distribution follow those of [Cerqueira et 


ah] (2006). For each extracted profile we take the radial 


^ocity of the peak, and calculate the difference (i.e., 
the radial velocity shift) for each pair of slits displaced 
symmetrically with respect to the jet axis: (3-5), ( 2-6) 
and (1-7). Each region (from I to IV; see Figure \T^ 
contributes to one radial velocity shift for each slit pair, 
and the mean value for a such radial velocity shift is 
plotted as a function of the distance from the jet axis 
in the right panel of the Figure (open circles; the er¬ 
ror bar gives the dispersion around the obtained mean 
value). In order to compare with the initially imposed ro¬ 
tational profile, we calculate for each radial distance the 
expected (at t = 0) radial velocity shift (for R > 0.2Rj): 
|Ai;rad| = 2 X 8'Cos0(i?j/i?), where the cos0 accounts for 
the inclination (red, solid curve in Figure p^. Although 
the match is reasonable for both inclinations (0 = 0° and 
0 = 15°), we can not reproduce the original profile for 
the case of 0 = 45°. The method applied here fails also 
to recover the original profile in the case where rotating 
and precession are combined, as in model M4. 
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6. CONCLUSIONS 


The PC A technique has been extensively applied in the 


astrophysical context (e.g., Steiner et ah 2009 


Ricci et ah 


2011 Malyshe^ 2012| ). PC A is used to identity quickly 


in large datasets patterns and correlations that other¬ 
wise would not be revealed. These patterns can then be 
interpreted in terms of arising from physically uncorre¬ 
lated phenomena present in the system, that thanks to 
PCA can be separated, identified and studied in detail. 
PC A offers also an efficient way to remove unwanted fea¬ 
tures from the data, as noise or instrumental fingerprints, 
and to compress and transmit only the relevant informa¬ 
tion present in a large dataset. We describe briefly the 
method in Section [2l The PCA method can be fruit¬ 
fully applied to the analysis of spectro-imaging data like 



observational data based on PCA tomograms is highly 
improved with the help of three-dimensional nu merical 
experiments. This has been shown, for examp l e, in|Heyer| 
fc SchloerbI (|1997|), Brunt, Heyer fc Mac Low (2009|) and 
ik & Blackman|(|2(Jlu|). in this paper we show 
)i the P 


Carrol, LranJ 
the potential applicability o: 


paper ^ 

CA to the study of HH 


jets by using the results of three-dimensional numerical 
simulations of rotating, processing jets. 

We apply the PCA to the synthetic images produced 
from the simulations, where we control the presence of 
known physical phenomena. We use the PCA to test 
if signatures of these phenomena can be isolated in the 
data, and produce a benchmark of PCA tomograms, with 
which PCA-processed real IFU observations of jets can 
be compared. In this work, we started generating from 
the simulations spectro-imaging datacubes around the 
wavelength of relevant emission lines for four different 
jet models. These are an intermittent jet model (Ml 
model), produced by imposing a periodic variation in 
velocity at the base of the jet, and models also including 
precession (M2 model), rotation (M3 model), or both 
rotation and precession (M4 model). For each of these 
models we applied the PCA decomposition to find the 
modes of variance with respect to the average intensity, 
and we discussed the tomograms corresponding to the 
modes of highest variance {k < 4). 

The tomograms, combined with their respective eigen- 
spectra, give important information on the processes oc¬ 
curring in the jet. The tomogram for k = 1 shows where 
the bulk of the emission varies the most, and it is in¬ 
dicative of the position of internal bow shocks generated 
by the imposed jet pulsation. The tomogram for /c = 2, 
in the precessing models, presents longitudinal gradients 
in intensity that we have identified with the precession. 
The slow variation of the intensity along the jet beam 
correlates with the jet orientation angle with respect to 
the observer, with abrupt changes in sign at the point 
where the jet inclination varies from directed toward the 
observer to away from the observer. The second tomo¬ 


gram therefore, can be used to determine if precession is 
present in the jet, and to quantify the precession angle. 
This feature is not present if precession is not included 
in the model. 

Tomograms 2 and 4 (for [O IJA6300 emission line), or 
3 and 4 (for [S IIJA6716 emission line) can be identified 
with the rotation of the jet around its axis in a rotating 
jet. We can use the reconstructed datacube to retrieve 
the initially imposed rotation profile. In a more general 
case, in which both precession and rotation are included, 
the signature for both is mixed in a given tomogram, 
and usually the signature is present in several different 
eigenvectors. For this reason, we are not able to retrieve 
the imposed rotation profile. 

Despite these promising results, it is worth to mention 
that there are some caveats to be kept in mind. In order 
to fully explore this method in real data, some further ef¬ 
fects should be carefully treated. The datacube should be 
as cleaned as possible and free of systematic effects. Slit 
positions should be placed along the jet beam, following 
the changes in the jet axis position due to changes in the 
jet propagation direction. The inclination of the jet with 
respect to the observer plays also an important role. It 
increases the number of relevant eigenvectors that must 
be taken into account. It also turns the analysis more 
complex, in the sense that several radial velocity gradi¬ 
ents appears in a given tomogram, making its interpre¬ 
tation more difficult. As the inclination angle increases, 
emission coming from the jet basis becomes more and 
more mixed with those from the internal working sur¬ 
face, which makes prohibitive the determination of the 
rotation profile. Also, the noise is expected to dominate 
the contribution to th e variance in the datase t at high 
order eigenvectors (see [Heyer fc Schloerb|1997 ), and real 
data will only give access to the first few eigenvectors. 

Nevertheless, the PCA remains a powerful diagnostic 
technique to analyze the structure seen in observation 
of jets in typical emission lines. We have shown that 
PCA can help disentangling precession and rotation in 
protostellar jets when the jet axis lies close to the plane 
of the sky. In a future work we will apply the PCA 
technique to the analysis of the HH 111 jet observed with 
the Gemini Integral Field Unit spectrograph. 


This paper was greatly improved thanks to the ref¬ 
eree, who has raised several and important questions 
during the review process. We are deeply indebted for 
his/her contribution to this work. We are thankful to T. 
Ricci and J. Steiner for enlightening discussions about 
the PCA technique and its applications to a datacube. 
AHC and MJV thank CNPq/CAPES for financial sup¬ 
port using the PROCAD project 552236/2011-0, and 
CAPES/CNPq Science without Borders program, under 
grants 2168/13-8 and 2565/13-7, respectively. We would 
like also to thanks Jerome Bouvier and Jonathan Eerreira 
for their kind hospitality during our stay in Grenoble 
(UJE/IPAG). EDG and JRI thank the DGAPA-PAPHT- 
UNAM (grants IA101413, IA103315). 


REFERENCES 

Anderson, J.M., Li, Z.-Y., Krasnopolsky, R. et al. 2003, ApJ, 590, Bacciotti, F., Ray, T.P., Mundt, R. et al. 2002, ApJ, 576, 222 

L107 Bacciotti, F., Ray, T.R, Eisldffel, J. et al. 2003, Ap&SS, 287, 3 

Bacciotti, F. 2004, Ap&SS, 293, 37 Bacciotti, F., Ray, T.R, Coffey, D. et al. 2004, Ap&SS, 292, 651 








































16 


Cerqueira et al. 


Blandford, R.D., & Payne, D.G. 1982, MNRAS, 199, 1982 
Brunt, C.M., Heyer, M.H. & Mac Low, M.-M. 2009, A&A, 504, 
883 

Cabrit, S., Pety, J., Pesenti, N., et al. 2006, A&A, 452, 897 
Cai, M. J., Shang, H., Lin, H.-H., & Shu, F. H. 2008, ApJ, 672, 
489 

Carrol, J.A., Frank, A. & Blackman, E.G. 2010, ApJ, 722, 145 
Casse, F., & Ferreira, J. 2000a, A&A, 353, 1115 
Casse, F., & Ferreira, J. 2000b, A&A, 361, 1178 
Chrysostomou,A., Lucas, P.W., & Hough, J.H. 2007, Nature, 450, 
71 

Cerqueira, A.H., Vasconcelos, M.J., Raga, A.C. et al. 2015, AJ, 
149, 98 

Cerqueira, A.H., Velazquez, P.F., Raga, A.C. et al. 2006, A&A, 
448, 2006 

Codella, C., Cabrit, S., Gueth, F., et al. 2007, A&A, 462, L53 
Coffey, D., Bacciotti, F., Woitas, J., et al. 2004, ApJ, 604, 758 
Coffey, D., Bacciotti, F., Ray, T.P., et al. 2007, ApJ, 663, 350 
Coffey, D., Bacciotti, F., & Podio, L. 2008, ApJ, 689, 1112 
Coffey, D., Bacciotti, F., Chrysostomou, A., et al. 2011, A&A, 
526, 40 

Coffey, D., Rigliaco, E., Bacciotti, F. et al. 2012, ApJ, 749, 139 
Coffey, D., Dougados, C., Cabrit, S. et al. 2015, ApJ, 804, 2 
Correia, S., Zinnecker, H., Ridgway, S.T., et al 2009, A&A, 505, 
673 

Davis, C.J., Berndsen, A., Smith, M.D., et al. 2000, MNRAS, 

314, 241 

Dougados, C., Cabrit, S., Lopez-Martin, L. et al. 2003, Ap&SS, 
287, 135 

Dougados, C., Cabrit, S., Ferreira, J. et al. 2004, Ap&SS, 292, 643 

Fendt, C. 2011, ApJ, 737, 43 

Ferreira, J. 1997, A&A, 319, 340 

Ferreira, J. 2008, New Astronomy Review, 52, 2008 

Ferreira, J., Dougados, C., & Cabrit, S. 2006, A&A, 453, 785 


Hartmman, L. 2008, in Accretion Process in Star Formation^ 
Cambridge University Press (2nd Edition) 

Heyer, M.H. & Schloerb, F.P. 1997, ApJ, 475, 173 
Kenyon, S.J., Dobrzycka, D., & Hartmann, L. 1994, AJ, 108, 1872 
Krasnopolsky, R., Li, Z.-Y., & Blandford, R.D. 2003, ApJ, 595, 
631 

Malyshev, D. 2012, arXiv: 1202.1034 

Melnikov, S.Yu, Eisloffel, J., Bacciotti, F., et al. 2009, A&A, 506, 
763 

Menezes, R. B., Steiner, J. E., & Ricci, T. V. 2014, ApJ, 796, 
LL13 

Mundt, R., Buehrke, T., Solf, J., Ray, T. P., & Raga, A. C. 1990, 
A&A, 232, 37 

Pesenti, N., Dougados, N., Cabrit, S. et al. 2004, A&A, 416, L9 
Pety, J., Gueth, F., Guilloteau, S., & Dutrey, A. 2006, A&A, 458, 
841 

Pudritz, R.E. 2004, Ap&SS, 292, 471 
Pudritz, R.E., & Norman, M.L., 1983, ApJ, 274, 677 
Pudritz, R.E., & Norman, M.L., 1986, ApJ, 301, 571 
Pudritz, R.E., Hardcastle, M.J., & Gabuzda, D.C., 2012, Space 
Science Reviews, 169, 27 

Raga, A. C., Navarro-Gonzalez, R., & Villagran-Muniz, M. 2000, 
RMxAA, 36, 67 

Ricci, T., Steiner, J., & Menezes, R.B. 2011, ApJ, 734, LIO 
Sauty, C., Cayatte, V., Lima, J.J.G., et al. 2012, ApJ, 759, LI 
Smith, M.D., & Rosen, A., 2007, MNRAS, 378, 691 
Solf, J., & Bohm, K.H. 1993, ApJ, 410, L31 
Spruit, H.C. 1996, in Evolutionary Proeesses in Binary Stars 
(R.A.M.J. Wijers et ah, eds.), pp. 249-286, Kluwer, Dordrecht 
Staff, J.E., Niebergal, B.P., Ouyed, R., et al. 2010, ApJ, 722, 1325 
Steiner, J., Menezes, R.B., Ricci, T. et al. 2009, MNRAS, 396, 788 
Testi, L., Bacciotti, F., Sargent, A.I. et al. 2002, A&A, 394, L31 
Vasconcelos, M.J., Cerqueira, A.H., Plana, H., Raga, A.C., & 
Morisset, C. 2005, AJ, 130, 1707 


APPENDIX 

EIGENVALUES 

In Table [2] we show the eigenvalues for the first ten eigenvectors, in terms of % of the variance, for the four models 
and the three inclinations. The first column displays the eigenvector E]^. The second column displays the inclination 
angle (j). The associated eigenvalues are show for each model (Ml to M4), for the computed emission lines [O I]A6300 
(columns 3 to 6) and [S II]A6716 (columns 7 to 8). 

The number of relevant eigenvectors, following the criteria of the scree test, is given in Table The scree test is a 
visual test used to obtain /Cmax, which is the maximum number of eigenvectors that we might consider in the analysis 
of the dataset. In this paper, we have defined a threshold for the eigenvalue, Athreshoid = 10“^ %, since we note that 
most of the eigenvalues will be intercepted by a line defined by such a value in a scree plot (that is, an eigenvalue 
versus eigenvector’s order plot; see Figure!^. This means that the rate of change in the eigenvalue as a function of the 
eigenvector’s order is almost constant (ana equal to zero, in our case) for /cmax < ^ < ^50, while it can varies a lot in 
the k < /Cmax range. With this in mind, we can look at Figurein order to find this point, that can easily recognized 
as the first point (from left to right) intercepted by the threshold eigenvalue function (black dashed line). 


THE EFFECT OF NOISE AND THE ENHANCEMENT FACTOR 


The effects of a noise background level in the data 

All figures in this paper have been obtained after a convolution procedure applied in each velocity channel map 
(VCM) in the real space (i, j, k) with a Gaussian profile, as explained in Section Compared with raw, non- 
convolved datacube, the net effect of the convolution process is to alter the value of the variance in the PC A treated 
datacubes. The first three columns of Table illustrate it, for a given model as an example (model M3, considering 
the [O I]A6300 emission line observed at an angle (j) = 0°). 

We have also degraded the data including in the simulation two effects: the presence of a random, background level 
of noise and considering also a broadened line profilj^ In the fourth column of the Table we show the eigenvalues 
for the first ten eigenvectors of such a model (again, considering the model M3, at 0 = 0° and the [O I]A6300 emission 
line). In Figure we present the tomograms/eigenspectra from k = 3 to k = 6 (from left to right, respectively). 
Different levels o^ignal to noise and inclination angles have been used: S/N ^ 10 and 0 = 0° (the first two rows of 
panels) and S/N ^ 5 and 0 = 45° (the last two rows of panels). 

The higher order eigenvectors that are dominated by noise can be disc arded in the reconstruction process, whic h 
constitutes on the other hand as a powerful filtering process for real data ( Steiner et al.||2009 Cerqueira et al.||2015 ). 


In the procedure to build the VCM’s we have artificially in¬ 
creased the FWHM of the Gaussian considered to spread the line 


intensity in radial velocity by a factor of ~ 2, in order to worsen 
the spectral resolution. 






PC A analysis of simulated jets 


17 


TABLE 2 



Eigenvalues, in 

% OF THE 

VARIANCE 

, OF THE 

FIRST TEN 

EIGENVECTORS. 


Eigenvector 

4> D 




Eigenvalue (%) 







[O I]A6300 



[S II]A6716 


Ek 


Ml 

M2 

M3 

M4 

Ml 

M2 

M3 

M4 


0 

99.7039 

90.3804 

88.3875 

75.5315 

99.7225 

97.4961 

92.6859 

84.1669 

El 

15 

89.8065 

79.8084 

71.6271 

53.2939 

97.2702 

96.4368 

85.6025 

79.3238 


45 

83.886 

76.2575 

62.3434 

49.7390 

94.6670 

90.2402 

80.3527 

70.0142 


0 

0.2597 

8.6249 

6.7645 

18.007 

0.2655 

1.8583 

3.9785 

12.7514 

E2 

15 

8.8301 

18.1922 

20.7925 

33.2961 

2.3652 

2.5940 

10.0743 

14.5136 


45 

10.7070 

16.6092 

22.8332 

33.4023 

3.6990 

6.3589 

14.4443 

18.6909 


0 

0.0339 

0.8386 

4.5483 

5.6351 

0.0108 

0.5996 

3.0637 

2.1859 

Es 

15 

0.8766 

1.3354 

5.1462 

8.9036 

0.2421 

0.7740 

3.5704 

4.6916 


45 

3.1367 

4.8132 

10.2354 

9.4788 

0.9499 

2.2272 

2.9921 

6.8108 


0 

0.0020 

0.1035 

0.2100 

0.6425 

0.0009 

0.0307 

0.2266 

0.7387 

Ea 

15 

0.2799 

0.3317 

1.5756 

3.3236 

0.0708 

0.0977 

0.4867 

1.0325 


45 

0.9187 

1.1429 

1.9413 

3.6313 

0.3297 

0.7967 

1.513 

2.8013 


0 

0.0002 

0.0424 

0.0822 

0.1289 

0(-5) 

0.0107 

0.0414 

0.1194 

E5 

15 

0.1300 

0.2048 

0.4929 

0.5452 

0.0298 

0.0575 

0.1233 

0.2162 


45 

0.7214 

0.4886 

1.2753 

1.4002 

0.1827 

0.1442 

0.2512 

0.8761 


0 

0.0002 

0.0068 

0.0031 

0.0326 

0(-5) 

0.0024 

0.0016 

0.0232 

Eg 

15 

0.0528 

0.0699 

0.2246 

0.3954 

0.0139 

0.0201 

0.0978 

0.1262 


45 

0.3256 

0.4470 

0.8301 

1.0521 

0.0921 

0.1339 

0.2133 

0.2717 


0 

0(-5)^ 

0.0020 

0.0029 

0.0166 

0(-5) 

0.0013 

0.0015 

0.0114 

Er 

15 

0.0133 

0.0364 

0.0755 

0.1059 

0.0041 

0.0109 

0.0209 

0.0468 


45 

0.1595 

0.0909 

0.2250 

0.6748 

0.0299 

0.0430 

0.1326 

0.2618 


0 

0(-6) 

0.0004 

0.0008 

0.0027 

0(-6) 

0.0003 

0.0003 

0.0015 

Eg 

15 

0.0062 

0.0120 

0.0411 

0.0811 

0.0021 

0.0052 

0.0152 

0.0297 


45 

0.0732 

0.0609 

0.1165 

0.2478 

0.0213 

0.0207 

0.0374 

0.1205 


0 

0(-7) 

0.0003 

0.0001 

0.0013 

0(-7) 

0.0001 


0.0007 

Eg 

15 

0.0019 

0.0045 

0.0128 

0.0318 

0.0007 

0.0017 

0.0050 

0.0101 


45 

0.0396 

0.0507 

0.0765 

0.1441 

0.0111 

0.0153 

0.0253 

0.0580 


0 

0(-7) 

0.0001 

-Of5) 

0.0003 

0(-7) 

0(-5) 

0(-5) 

0.0002 

Eio 

15 

0.0014 

0.0017 

0.0061 

0.0106 

0.0003 

0.0009 

0.0019 

0.0037 


45 

0.0163 

0.0220 

0.0536 

0.0942 

0.0098 

0.0111 

0.0140 

0.0407 


^The order of magnitude will be indicated instead of eigenvalue if it is less than 10 


In this particular case, 0(-5) implies that A 


10 -^%. 


TABLE 3 

The number of eigenvector that should be retained in the data considering the scree test. 


Model 

[O I]A6300 

^scree 

[O I]A6363 

test 

[S II]A6716 

[S II]A6731 



4 

4 

4 

3 

0 

Ml 

7 

7 

7 

4 

15 


11 

11 

11 

9 

45 


6 

6 

5 

4 

0 

M2 

8 

8 

7 

6 

15 


11 

11 

11 

9 

45 


6 

6 

6 

6 

0 

M3 

9 

9 

9 

9 

15 


13 

13 

12 

9 

45 


8 

8 

7 

7 

0 

M4 

11 

11 

9 

8 

15 


14 

14 

13 

11 

45 


The enhancement factor 

The tomograms represent an image in a new system of uncorrelated coordinates. We can use the tomograms, 
together with their eigenspectra, to identify different physical properties present in the data. We can then reconstruct 
the datacube, omitting some undesirable property and putting in evidence others. In particular, we want to identify 
and discriminate from the tomograms calculated so far, the presence of precession and rotation. In the Section we 
have shown several tomograms and eigenspectra. 

Once identified, a given property “A” can be emphasized using the equation (§: 


I'p^{A)=T0kikA)-[ExkikA)f, 


(Bl) 


to reconstruct the datacube. In this process, we keep the eigenvectors already identifie d and associated with the 
property “A”, in the characteristic matrix. This can be done following the prescription in Steiner et al. (2009): we 
can multiply the columns of the characteristic matrix by a factor, T/^, which can be 0 or 1. It acts suppressing a given 
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eigenvector (or Tk = 0) and keeping others (F/. = 1) in the characteristic matrix. The end product is a matrix that 
has the desired property only (feature enhancement). 

In the Section we have shown the tomograms associated with a specific eigenvector, and extracted information 
about the properties they contain. In Section the reconstruction process has always been guided by these findings, 
where the tomograms and eigenspectra have been analyzed together in the search for the precession and/or rotation 
signatures. 


TABLE 4 

Eigenvalues, in % of the variance, for the first 10 eigenvectors for model M3 (rotating model), for the raw, convolved 

AND NOISY DATA. THE [O I]A6300 LINE HAS BEEN CONSIDERED FOR (p = 0°. 


Eigenvector 

Ek 

raw data 

Eigenvalue (%) 

+ convolution 

+ noise 

El 

99.190 

88.3875 

99.182 

E 2 

0.7523 

6.7645 

0.481 

E 3 

0.0572 

4.5483 

0.057 

E 4 

0.0033 

0.2100 

0.009 

E 5 

0(-5) 

0.0822 

0.009 

Eq 

0(-6) 

0.0031 

0.009 

Er 

0(-7) 

0.0029 

0.008 

Eg 

0(-8) 

0.0008 

0.008 

Eq 

0(-9) 

0.0001 

0.008 

Eiq 

O(-IO) 

0(-5) 

0.008 
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Fig. 19.— Sequence from the third (left) to the sixth (right) tomogram and eigenspectra for M3 model and [O l]A6300 emission line with 
a background noise. The two topmost panels are for S/N ~ 10 and 0 = 0° and the bottommost ones are for S/N ~ 5 and 0 = 45°. The 
contribution to the variance is mainly dominated by noise fluctuation in high order eigenvectors. 





























































